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Effects of spatially varying interfacial parameters on the propagation of surface waves are studied. 
These variations can arise from inhomogeneities in coverage of surface active substances such as 
amphiphillic molecules at the fluid/gas interface. Such variations often occur in phase coexistence 
regions of Langmuir monolayers. Wave scattering from these surface inhomogeneities are calculated 
in the limit of small variations in the surface parameters by using the asymptotic form of surface 
Green's functions in the first order Born approximation. When viscosity and variations in surface 
elastic moduli become important, modes other than transverse capillary waves can change the 
characteristics of propagation. Scattering among these modes provides a mechanism for surface 
wave attenuation in addition to viscous damping on a homogeneous surfactant covered interface. 
Experimental detection of waves attenuation and scattering is also discussed. 
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PACS numbers: 68.35.Ja, 47.35. +i, 92.10.H, 68.10.-m 



I. INTRODUCTION 



Investigation of interfacial phenomena has attracted 
much attention, especially if interfacial properties are 
modified by a surface active material. Many studies 
involve amphiphillic molecules at the fluid/air interface 
with monolayer or less coverage. Small amounts of mate- 
rial at an interface can drastically alter wave propagation 
and also lead to complex phase behavior. Recent experi- 
ments with fluorescence and Brewster angle microscopies 
have sought to elucidate the structure of monolayer films 
of fatty acids, cholesterol, and phospholipids. H A myr- 
iad of phase transitions of the surface film |are directly 
observable as well as theoretically predicted.a 

These studies often reveal high surfactant concentra- 
tion domains of various shapes and sizes in the coexis- 
tence regions. Figures 1(a) and 1(b) show Brewster an- 
gle microscope images of circular and striped domains of 
triglycerides in the Liquid Condensed/Liquid Expanded 
and Liquid Expanded/Gaseous "coexistence regions" re- 
spectively. The figures span about Imm. Lighter re- 
gions correspond to regions of higher surfactant concen- 
trations. This type of coverage is often associated with a 
phase coexistence region not unlike that of the gas-liquid 
phase transition. These first-order like transitions can be 
described approximately by a thermodynamic Maxwell 
construction across mechanically unstable regions of the 
equation of state. External electric fields can also be 
used to generate local interfacial surfactant concentra- 
tion gradients. H'Q If nonlocal interactions such as dipole 
forces are important, regular modulated fihases appear 
in the usual two-phase coexistence region.B 

Experiments have also investigated changes in capil- 
lary wave propagatiorLcaused by insoluble surfactants at 
an air- water surface. Q'El Wave damping experiments on 
film covered air/water interfaces show a remarkable in- 
crease in the damping rate above that caused by viscous 



bulk fluid dissipation. Even a single monolayer of sur- 
face active substances, such as amphiphillic fatty acids, 
increases the damping rate by as much as five fold by 
altering the flow trajectories beneath the surface. In ad- 
dition, at monolayer concentrations, the damping coeffi- 
cients appear to be sensitive to nonuniform surface cov- 
erages. At concentrations where an apparent coexistence 
of surfactant phases occurs, the damping rate peaked at 
five or six times that due to viscous dissipation. Q 

Attempts to explain the propagation of water waves 
covered with surfactant have considered the homoge- 
neous mechanical properties, of fioatuig filmsEI. Althouffih, 
the theories of Dorrenstein,Ll Levichjj, and LuccasseroH 
all predict a maximum in the damping factor, they are 
incomplete since they rely on surface parameters which 
are not well defined in a nonuniform surface. Though a 
peak in damping rate is experimentally observed, these 
systems are almost always inhomogeneous. As argued in 
Refs. (11) and (12), wave propagation should depend on 
the heterogeneous nature of the interface. The details of 
this surface inhomogeneity depend on how the fllm was 
applied as well as the phase diagrams of these surface 
materials. 

Wave scattering off of surface heterogeneities can serve 
as an additional damping or attenuation mechanism. 
Much as light is attenuated in a turbid medium, sur- 
face waves can be scattered by inhomogeneities like those 
in Fig. 1. Gravity wave scattering from a (discontinu- 
ity in the trough depth has been calculatedll3 and ex- 
tended to multiple scattering-and wave localization in a 
random bottom by Belzons.llil Single scattering of grav- 
ity waves o£u3bstacles or depth discontinuities has been 
well studiecO and has been applied to multiple scattering 
of gravitA^ waves off obstacles protruding from the water 
surfacdlj. Work done in this area has usually been in 
the context of oceanography and has focussed on gravity 
waves propagating in the presence spatial variations of 
the boundaries, particularly the finite depth bottom. 



We will consider capillary waves scattering off varia- 
tions in parameters such as surface elastic constants and 
surface tension. An exact solution of water wave scatter- 
ing from discontinuities in boundary [Condition parame- 
ters can be found for very special casesll3; in this paper we 
will use approximation methods which are valid for small 
variations in surface tension and interfacial modulus. 

In the next Section, we present a model for an inviscid, 
infinite depth substrate fluid with a spatially varying sur- 
face tension. Regions of sharply different surface tensions 
are stabilized by line tensions at their boundaries. This 
line tension could be supplied, for example, by a waxed 
circular thread floating at the interface and separating 
a surfactant rich monolayer from a clean fluid interface. 
This model problem leads naturally to a Green's func- 
tion describing the response due to a point normal force 
acting on a uniform surface. 

In Section III the asymptotic limit of this Green's func- 
tion is used in a Born approximation to describe the scat- 
tering of potential flow by a circular domain of different 
surface tension. The scattering function is used to cal- 
culate the attenuation of a plane capillary wave via a 
scattering cross section. 

The effects of bulk viscosity and surface viscoelasticity 
are included in Section IV. Variations in surface com- 
pressibility are the primary source of scattering in this 
context since the surface tensions of coexisting island 
phases like those in Fig. 1 are nearly equal. Here the 
tangential stresses mix both potential and vorticity fields. 
The lowest order Green's function matrix is derived. The 
analysis shows the existence of at least three distinct sur- 
face modes: transverse, "in-plane" longitudinal, and "in- 
plane" shear. 

Scattering via the Born approximation for a viscoelas- 
tic surface is calculated in Section V using the Green's 
functions of Section IV. The dissipation of energy into the 
various channels is also considered. Wave amplitude and 
energy attenuation due to these processes is discussed 
in Section VI, where experimental implications are also 
considered. 



small, vanishing as the reciprocal island dimension for 
the configurations shown in Fig. 1. Allowing for fixed 
surface tension variations is nevertheless instructive and 
might even be realizable experimentally over a reason- 
able time period by confining one phase of an insoluble 
monolayer inside a flexible barrier of negligible stiffness. 
The mechanical stress balance at the interface in 
the psesence of a spatially varying surface tension a(x) 
reads ,tZl 

n • (n+ - n_) = a{r} (^-L + -L^ fi - V. a{r), (2.1) 

where n is the unit vector normal to the interface, V^ is 
a tangential gradient on the interface, and Ri and i?2 are 
the local principle radii of curvature of the interface, rib- 
and n_ denote the stress tensor above and below the 
interface respectively. On an inhomogeneous surface, 11 
and a may depend on the surface coordinates; the projec- 
tion of these coordinates onto the x — y plane is denoted 
by r. For ideal fluids, the stress tensor is proportional to 
the pressure. In the Monge representation, the linearized 
normal component of Eq. (ET[) is. 



7r(z =: 0) = cr(f)Vl77(f), 



(2.2) 



where 77 is the vertical displacement of the interface from 
its equilibrium positio n at z = 0, and V^ is a gradient in 
the X — y plane. Eq. ( ^.2| ) is the usual dynamic bound- 
ary condition which determines the dispersion relation 
for surface waves. Here, tt is the jump in dynamic pres- 
sure across the interface. We assume that the tangential 
stresses are balanced by surface forces such as line ten- 
sion. 

For an ideal fluid, the motion will be irrotational and 
the velocity field can be expressed via a potential. 



V'l = di4>{f, z) 



(2.3) 



satisfying V^(/)(r, z) = in the bulk fluid. At the surface 
we have, 



II. MODEL FOR SURFACE TENSION 
VARIATION ON AN INVISCID FLUID 



In this section we investigate capillary wave propaga- 
tion in the presence of surface tension inhomogeneities 
with an inviscid Newtonian fluid occupying the half space 
z < 0. This first model is somewhat artificial because if 
one hypothesizes a surface tension discontinuity due to 
two distinct surfactant phases at an interface, one phase 
would immediately grow at the expense of the other to 
reduce the free energy. With short ranged interactions 
only, the surface tensions of two phases coexisting in 
equilibrium are necessarily equal. The surface tension 
changes across modulated phases with large scale struc- 
tures determined by dipole-dipole interactions are also 



dtTiif)~d,cl){r,z = 0), (2.4) 

and the Navier-Stokes equations lead to, 

7r(f, 2) = -pdt(f>{f, z) ~ g-q{r), (2.5) 

where 7r(r, z) is the fiuid pressure, p is the fiuid density 

and g the gravitational constant. Henceforth, gravity 

will be neglected [g — Q). The equation for the velocity 
potential at the surface reads, 



[-dl + (ao + 5<j{r))yld,] c^{r, z) 



z=0 







(2.6) 



We now introduce a harmonic time dependence w to 
all quantities and set yO = 1. Upon adding a unit 
z— component delta function force to Eq. (2.2), we are 
led to define a surface Green's function which solves 



[^2 + [ao + Sa{f))Vld,] G{f, 0) = -^iuj5{r). (2.7) 

The deviation of the spatially varying surface tension 
from the constant value ctq is denoted by 6a (f) . Equa- 
tion (2/7) can be solved for Sa = Ohy expanding Go in a 
general Fourier integral, 



Go{f,z,u) 



(2^)2 



Go(fc)e*^■*'+l'=l^ 



(2.8) 



The solution to (2^) can then be formally written as, 



G(fc) = Go(fc) + Go(fc) 






q^\q\5a{k-q)G{q) 



with 



aoGo(fc) = (fc2|fc|_p3)-i 



(2.9) 



(2.10) 



where p = (w^/cro)^ is real and positive. This trans- 
verse mode described by Go{k), i.e., capillary waves prop- 
agate at fc = zLp. 



III. SCATTERING FROM A SINGLE CIRCULAR 
DOMAIN 

In this section, the scattering of a plane inviscid capil- 
lary wave from a single circular domain with a sharp sur- 
face tension discontinuity is calculated. The surface ten- 
sion inside a circular domain of radius a is cti = ctq + ^f- 
A perimeter whose only effect is to confine the domain 
and maintain the discontinuity is assumed. A schematic 
of this scattering problem is depicted in Figure 2. The 
wedge at the left is typically -an electrocapillary wave- 
maker, excited at frequency ujIlj 

Two approximations will be made. First, a transverse 
capillary wave impinging on a domain will change its 
perimeter; we neglect line tension energy of the confin- 
ing loop in comparison to the areal stretching energy of 
the surface. Second, we will assume a frozen variation in 
surface tension due to a fixed domain. This last simplifi- 
cation requires weak scattering, i.e., ka <ti 1, where k is 
the wavevector and a is the domain radius. 

For a single domain, a two dimensional analog of the 
Young-Laplace equation can be defined for circular do- 
main of radius a with a sharp surface tension discontinu- 
ity. 



CTO — CTl 



2. 

2a 



(3.1) 



where ctq and CTi are the surface tensions inside and out- 
side the domain respectively and 7 is the line tension 
actin g through the circular domain boundary. Eq uation 
(3.1) is actually the tangential component of (2.1) when 
bulk viscosity is neglected. For small Sa, the scattering 
is obtained in the first iteration or Born approximation. 



In the inviscid case, the elastic scattering mechanism de- 
scribed here will be the only process resulting in attenua- 
tion of an incident plane wave in the midst of a collection 
of scattering centers. 



By using Eq. ( 2.10 ) , we can find the exact surface 
Green's function, Go{f), for a homogeneous interface, 



Goir) = ^ 



d^k 



-ik-r 



aoJ (2^)2 fc2|fc|_p3_ 



i^"™'-) 



TT^aop Jq 



yd +pe 



dy (3.2) 



In the above integrands, p+ = p + ie represents a small 
damping which removes the ambiguity in the integrals 
for outgoing Green's functions. Since p is real, ( p^ ) is 
dominated by the first term at large distances {pr 3> 1) 
and the asymptotic form becomes. 



Goir) 



^ipr+37Ti/4 



3aoP ^/27^pr 



(3.3) 



With the help of the Green's function, Go{r), the solution 
of Eq. (2.7) can be written in integral form as, 



(r, 0) = e'P" + - / d'r' ,5a(r') Go(|r - f'|)V "9.0(r', 0) 

(3.4) 



The first term of (3.4) represents an incoming plane 



wave with unit magnitude at the origin of the scat- 
terer generated by a z component force at a; ^ —00 
such that the power law evanescent wa ves represented 
by the second term in the last line of (3.2) have long 
died off. The integral term is the scattered wave due 
to the "potential" Sa{r'). In the lowest order approxi- 
mation for small (Scr(r'), 0(r') in the integral is replaced 
with <j>incir') = e^P^ . In the far field, the potential takes 
the form. 



Hr) 



-Jr 



(3.5) 



where the s catt ering factor is the asymptotic limit of the 
integral in (^.4[), 

Validity of this approximation is ensured when the scat- 
tered wave just outside the disc is small compared to the 
incident plane wave, which leads to the requirement. 



{paf 



<1. 



(3.7) 



In this asymptotic limit, the total scattering length is. 



E, = 



271- 



d0\fid)\' 






(3.8) 



Upon using the capillary wave dispersion formula, ui = 
^/oqp^, we see that Eq. ( |3.8| ) imphes that E^ ~ w^, in 
contrast to the famous RayleigJMesult S^ ^ tj^ for light 
scattering from the atmosphereta. 

Nonuniform coverage of air/liquid interfaces however, 
usually results in a collection of circular scatterers. Scat- 
tering from an uncorrelated collection of identical do- 
mains (assumed for simplicity here to be of identical size) 
reduces the forward coherent intensity according to, 



|(0(x))|2oce-"^-'(P)" 



(3.9) 



where n is the numb er o f scattering domains per area. 
Beer's law, equation (3^), applies in the single scatter- 
ing regime and is valid when the typical inter-domain 
distance I 3> 27r/p 3> a; our calculation for E^ ^ a 
only describes weak attenuation accurately. An alterna- 
tive calculation which averages over surface disorder and 
leads to the same result is given in Appendix B. 



IV. MODEL WITH BULK VISCOSITY AND 
SURFACE ELASTICITY 

In actual experimental systems such as the one de- 
picted in Figure 2, surface layers have additional param- 
eters which are spatially varying. For example, the sur- 
face compressibility will be different inside and outside 
the islands even if the surface tension is uniform. The 
restoring forces proportional to these surface parameters 
are coupled to the bulk viscous forces at the interface. 
Surface elasticity is an important factor in the propaga- 
tion and attenuation of surface wavea3 and can arise from 
surfactant molecules at the interface. Fluid motion in the 
interfacial plane tends to locally distort the film which re- 
sists due to its surface compressive or shear elasticity and 
interfacial 2D viscosity. 



The normal component of (2.1) now reads 



(do + (5CT(r))V^U2 — iujn{r, z) + 2iu)vdzVz 



(4.1) 



where a^ is the uniform background surface tension, 
5a {r) is the spatial variation in surface tension consid- 
ered earlier, v is the bulk fluid viscosity, and fzi'r) is 
proportional to an external pressure in the ^-direction. 

Upon considering the intrinsic stresses of the interface, 
we find that the tangential forces are balanced when, 



V_icr(r) + v_i-n 



(2) 



y{dzvi_ + V^Vz) 



/±(r), 



(4.2) 



where /j^ (r) is proportional to an in-plane stress ap plied 
at the interface. The first term on the left side of (4.2) 



arises from gradients in surface tension and acts like a 
pressure term in the stress tensor of an ideal two dimen- 
sional liquid. The dissipative part of the 2D fiuid stress 
tensor is given by 11 ^"^K In the remainder of the paper, 
we assume V_l • n ^^^ — 0; all of the phenomena can be 
adequately accounted for by Vj_cr(7^). Modifications due 
to dissipative and elastic contributions to 11 '^^ is consid- 
ered in Appendix A. 

Fluid motion at the surface convects surface molecules 
and changes the surfactant concentration r(r, t) locally, 
so we set 



^Mr) = Jv^r(f). 



(4.3) 



For an insoluble niaterial, conservation of surfactant at 
the interface gives,! 



dtn^ + V^-{v^T{r)) 



'iuj5T{r) + TQV^-v^ =0, 

(4.4) 



with r(r) = Fq + 5T{r)^ where Fq is the equilibrium 
surfactant concentration inside or outside the domain. 
Equations ( |4.3| ) and (4.4) can be combined to express 
the surface tension gradient in phase j in terms of the 
in-plane velocity field. 



V_LCT(f)~ -Bj\I^{V^-vi_) 



where the bulk modulus is defined by 
„ da 



B, 



'dV 



r=r,- 



(4.5) 



(4.6) 



with j = 0, 1 representing B outside or inside a do mai n 
for example. We can now solve equations (4.1) and (p3). 
Upon decomposing the velocity into its potential and vor- 
ticity components, 



Vi = di(t) + ei-jkOji^k 



(4.7) 



we choose a gauge where tpz 
(the vorticity is fJ = V x i/ = 



V0(r) 



0, and define j3i = etji^j 
-z X V^/3) so that. 



dzfi 



(4.8) 



where the vector /3 lies in the x — y plane. The fields 
and /3 satisfy. 







(4.9) 



and 



dt(3 = iyV^f3 (4.10) 

respectively, which are solved by the expansions. 



_fk_ 

(2^ 



■(j){k)e 



ik-r-\-\k\z 



and 



/3(r,z) 



(2^ 



m 



ik-r-\-lz 



(4.11) 



(4.12) 



gential stress balances at the surface z ~ after Fourier 
transforming in the in-plane r coordinate take the form, 



Mq)lP{q) = F{q) 



_fk_ 



M{k,q)-ip{k) (4.13) 



where I = Xjk"^ — iuj /v. The positive roots for / and |fc| 
are taken to ensure th at (p a n d (3 v an ish as z — > — cx j. 
Finally, upon using ( pi] ), ( |til ), ( ^ , (fj) in (p]), 



and taking a time derivative of (4.2), the normal and tan- 



M{k,q) = - 

UJ 



where F{q) is the g*'' Fourier component of {fz, fx, fy)- 
In ( 4.13 ), tpir) = {(f>, /3x, f3y) and k and (fare wavevectors 
in the x — y plane. Spatial variations in the surface ten- 
sion (j{r) and interfacial bulk modulus B are included in 
the matrix 



(4.14) 



6(T{q- 


-k)P\k\ Sa{q-k)k^kx Sa{q-k)k^ky 


5B{q- 


- k)iPkx SB{q- k)kll SB{q- k)kxky 


5B{q^ 


- k)ik^ky 5B{q~ k)kxkyl 5B{q- k)kll 



When the film is homogeneous, M = 0. 

The coefficients a and B in L(g) are constants, e.g. the values uo and Bq outside an isolated circular domain, or 
mean values appropriate for a heterogeneous collection of many domains. In equilibrium, for large domains, the surface 
tension contrast will be small as discussed above. Hence, in the following analysis (unlike that for the nonviscous 
fluid), we make the simplifying assumption that the surface tension is constant in monolayers in coexistence and near 
equilibrium. Scattering will be solely due to the contrast in B. The compressional modulus for phase j = (0, 1), Bj, 
can be measured by the slopes of the isotherms just outside the coexistence plateau as shown by the dotted lines in 
Fig. 3. These values will be different in the two phases even for monolayers in equilibrium. 
The matrix L((j) corresponding to an interface with B = Bq is, 



Mq) 



J — 


^\q\ + 2tuq^ 


iqx 

iqy 


■f^?' + 2H9I 

['^''q^ + 2u\q\] 



Qx 



2wl 



cpq 



qy 






2ivl 



<^oq 



[^^'yl + HI' + ll)] 



(4.15) 



Upon inverting (4.13), the equation for xp takes the form 



TPiq) = Giq)F{cf) + G{q) 



d'k 

(2^)2 



M{k,q)rp{k), (4.16) 



with a matrix Green's function G{q) = L~-'-(g). For fixed 
UJ, the poles of the Green's function give the wavevectors 
of the possible modes. For the region of parameter space 
most accessible to experiments, there are two poles corre- 
sponding to outgoing capillary waves. Appendix A shows 
that if Vj^- n '^^^ ^ 0, additional modes corresponding to 
in-plane motion also exist. 

When viscosity is small, the "transverse" capillary 
wave pole is near the inviscid pole; it is shifted from k = p 
by a small real part and develops a small imaginary part 



due to viscous damping. When x = i' (cj/ctq) <C 1, 



The functions,EHij 




3 V Co / j/^ — y v2a; + X 



(4.17) 



(4.18) 



and, 



A" 



^/^ y2 ^y^ _ 2yxV2x + Ax^ 



y\/2x -I- X 



(4.19) 



are plotted in Figure 4(a) as a function oi y = Bo/ao- 
For LU <C B"^ /v^ the longitudinal pole is approximately. 



,1/4^3/4 

Vb^ 



,i^/8 



(4.20) 



The relative dissipation between the transverse and lon- 
gitudinal modes, Im{kt\ / I'm{ki\ , are plotted in Figure 
4(b). The "shear" mode found in the Appendix A can 
usually be omitted since it has a large imaginary part and 
is highly damped. The asymptotic form of the Green's 
function including the transverse and longitudinal modes 
is explicitly displayed in Appendix C. 



V. FIRST ORDER SCATTERING FROM 
VISCOELASTIC DOMAINS 

To consider the effects of surface elasticity, the inclu- 
sion of a nonzero bulk viscosity is necessary to satisfy 
the tangential force balance equation. Thus, there is in- 
elastic wave attenuation even without surface inhomo- 
geneities and elastic wave scattering. In this Section, we 



calculate the scattering of surface waves from domains of 
different modulus B only. The variation in the modulus, 
6B = Bi — Bq, oi domains near equilibrium can be es- 
timated from the isotherm construction in Fig. 3. Since 
the periods of the surface waves are probably smaller 
than the molecular diffusion times, the relevant quantity 
is the isoentropic modulus rather than the isothermal 
modulus found from Fig. 3. We assume equality of sur- 
face tensions inside and outside the domains. 

Incident waves with various amounts of transverse and 
longitudinal (and possibly shear) mode components can 
be generated with surface forces at various angles from 
the normal. The incoming wave we consider will be that 
generated by a force which acts solely in the z direction 
at the surface. The relative amounts of transverse and 
longitudinal components is fixed by the coefficients of the 
relevant Green's functions given in Appendix C. 

In the large r limit, Eq. (4.16) yields the form. 



■0(r) ~ xp„{x) 



E 

i.b=t,l 



9a.,m- 



r,ika.r 



(5.1) 



where iPq{x) is a plane normalized combination of trans- 
verse and longitudinal waves generated by a force per 
length fz with infinite extent in the y direction. The wave 
induced by the wavemaker will initially have transverse 
and longitudinal modes with both and /3x components. 
By integrating the asymptotic forms of the Green's func- 
tion from Appendix C with respect to the direction along 
the wavemaker blade, we can determine ipQ to leading or- 
der in ly. The vertical interfacial displacement generated 
by a plane wavemaker. 



ri(x) 



^ikfX 



3croP 



e^"/^^^/% .fe,. 



2b;]/'wi/4 



(5.2) 



determines the nature of the waves incident on the scat- 
tering domains. If Im{kt} <C Im{ki}, incident waves are 
mainly transverse capillary waves; if Im{kt} 3> Im{ki}, 
the incident waves are predominately longitudinal waves. 
The ratio of the imaginary parts of the transverse and 
longitudinal modes are plotted in Fig. 4(b) for various 
values of B. 



The form factor in Eq. (5.1) describes a complicated 
scattering process. The 4> and (3 components of xp scatter 
among themselves; as suggested by the indexes a,b = t,l 
on the scattering amplitudes ga^b{Q)- The components of 
ga.b{0) governing the scattering of (j) has four terms de- 
scribing scattering between transverse and longitudinal 
modes. To lowest order in iov^ /a^ and ojv^ / B^ these 
are. 



9t,m 



Aa' 



■nvp 



,7/4 



1/4 



t/4 



6B\ Ji{2apsinT 



Bo J 2apsin^ 



(5.3) 



5m(^) 



e5-/8a2^p 



2 ^^i7/8i'2Bo~ao)iy^/^ fSB 



Bq CTq 



3/2^1/8 \Bo 

Ji (qa) 



qa 



cose, (5.4) 



m 



i{e)^e'^'^y/2^i{pa)- 



3cro(2Bo-cro)i^^ (5B 



ABl \Bo 



Ji(2afc;sm|) 
2akisin^ 



cose (5.5) 



and. 



,., 2 9/4 f^ To^'i^'/' fSB\ J,[qa) 
9i,t{^) ^ a p-'' \ t:^-^ — 1 -fT ) —— ^os0. (5.6) 



2 B^ \BoJ qa 



In the above equations, q = ^Jk'f A- k'f — 2ktkicose. Each 
of these scattering functions correspond to a particular 
intermode scattering. For example, gt,i{0) represents the 
scattering of an impinging transverse wave into a longi- 
tudinal wave. The two limiting cases depicted in Figures 
5(a) and 5(b) will be of interest. 

In Fig. 5(a), the longitudinal mode is more dissipa- 
tive than the transverse mode. For modes described by 
Eqns. (pT7|), (|t1), (pT9|), and (Ooh, this first case 



corresponds to Im{ki\ ^ l~^ ^ A". In this case, we 
may assume that only transverse capillary waves reach 
the scattering centers. Some of the transverse amplitude, 
however, is scattered into the longitudinal mode which is 
then quickly viscously dissipated leading to enhanced in- 
elastic attenuation of the velocity potential. For an areal 
density of identical island scatterers n, the transmitted 
potential decays according to £3 

Tl 

(|0(a;)|) -- exp[--{Y.t,t + ^t,i)x - ax] (5.7) 

where the scattering and absorption lengths, Ej^j and Tit,i 
are given by angular integrals over gt,t{e) and gt,i{Q), 



A{pafTi'^v fSB 
B^ 



(5.8) 



and, 



^t,i 



{pa)\^jy^/V/^ 



i2Bo-aof{^-£] . (5.9) 



To make these and subsequent expressions dimensionally 
meaningful, let ao, Bq —>■ ao/p, Bq/ p where p = 1 g/cm^ 
for the interface between air and water. Since uj ~ p^''^ 
for capillary waves, these scattering lengths depend on 
frequency and viscosity according to Y^t,t ~ uuf'^, and 
St./ ~ t/^/2^17/6 ]\jote that E^^j describes an elastic scat- 
tering process, which could be detected by placing a wave 
detector at an angle to the incident wavevector, and mea- 
suring the scattering amplitude gt,ti(^) directly. J^t,i acts 



as an enhanced viscous dissipation under the conditions 
described above. In our approximations, Ej^j ^ St /; the 
major scattering process is elastic. Only when (Tq/Bq is 
very large can T,t,i be comparable to Sf^t. The damping 
from vis cous dissipation is described by a = A" displayed 



in Eq. (|J_ 



In Fig. 5 (b), the opposite situation obtains, A" > 
Im{ki\. The predominant mechanism is scattering of 
longitudinal waves. Here, the relevant scattering lengths 
are. 



^l^l ~ 9^2(p^)4 



^i«/y/8^17/4(2B„_^„)2 



32B, 



13/2 



5B_ 
Bo 



and, 



^i,t 



TT^ipaf 



5/2 3 



8Bq \ Bo 



(5.10) 



(5.11) 



which behave as E;^; ^ t/^'^'/^a;'*^/^^ and E;^t ~ (z/w)'^. 
The approximations oji'^ /Bq <C 1 and uiv^ ja^ ^ 1 give 
Ej^t ^ E/_;; in this case, the predominant scattering 
is that of longitudinal waves scattering into transverse 
waves which are then quickly dissipated by the bulk vis- 
cosity. Again, only when uojBo ^ 1 is the elastic lon- 
gitudinal — *■ longitudinal scattering comparable. Here, 
simple viscous attenuation of longitudinal waves behaves 
as, 



a = Im{ki} 



\/2W2 Z.l/4^3/4 



(5.12) 



VI. IMPLICATIONS FOR EXPERIMENTS 

Surface wave damping is expected to have experimen- 
tal signatures with the techniques used in recent studies. 
Induced surface waves can be probed by light scatter- 
ing and detecting, the beam deflection or measuring the 
power spectrumJ2l Alternatively, direct mechanical mea- 
surements due to the force from impinging waves are also 
used. 

Our calculations for one scatterer can be tested by 
scattering an induced wave off a single isolated domain. 
The asymptotic radial distribution of scattered poten- 
tial should follow those given in equations (3^) and (5.1) 
for the elastic (inviscid) and inelastic (viscous) cases re- 
spectively. Because a single surfactant domain is rarely 
seen, an isolated domain may have to be constructed by 
placing a thin confining thread on the surface of a ho- 
mogeneous monolayer and aspirating the material out- 
side of this loop. Provided the solubility into the bulk 
liquid is small, this arrangement, although metastable, 
may be quite long-lived. In practice, the scattering from 



the thread loop alone may have to be considered. Elec- 
tric fields have also been used to manipulate surfactant 
concentrations and orientations at the interface.0a Sin- 
gle circular domains of surfactant rich or depleted phases 
have been produced with this mcthod.Q 

Damping of induced height fluctuations propagating 
through an array of scattering domains can be directly 
measured by probing the amplitude as a function of dis- 
tance from the source. Experimentally, the nucleation of 
surfactants results in domains of various sizes and shapes. 
Our analysis describes a monodisperse collection of scat- 
terers but can be extended to include polydispersity by 
averaging over domain size distributions. 

Using Eqs. ( |3.8D and ( |3.9| ) as guides, an estimate of the 
scattering effect can made. For a membrane which has 
only a surface tension variation, scattering attenuation 
above that of viscous damping will be noticeable when. 



> 



a{pa)^ao 



(6.1) 



The value Sa is difficult to measure directly. The surface 
tension difference inside and outside a domain of radius a 



is related through Eq. (3.1); however, these quantities are 
difficult to measure but would be expected to be small for 
systems near equilibrium. The plateaus in the isotherms 
are rarely perfectly flat due to long-ranged dipolar in- 
teractions between surface molecules and the slom-relax- 
ation to equilibrium. An experimental estimatecEl gives 
line tensions on the order of 7 ~ 10~^ erg/cm. For do- 
mains of radius r^ 0.1 mm, this tension gives Sa ~ 10"'^ 
dyne/cm. For a fluid with cto — 70 dyne/cm driven at 
a frequency w = 27r x 100 s~ ^ and covered with ^^ 10'^ 
scatterers per cm^, Eq. ( |6.1[ ) requires unphysical bulk 
viscosities v < 10^^"^ cm /s before attenuation due to 
scattering will be noticeable. 

When viscosity is explic itly included, the transverse 
wave attenuation given by (5.7) includes a damping term 
—A"x which has a maximum (Fig. 4(a)) as predicted by 
homogeneous theories of Dorrenstein and others.LrQ For 
the coexistence of fatty acids in the low coverage regime, 
the modulus B is expected to be much smaller that the 
surface tension; thus is not clear whether the maximum 
seen in experiments can be described by this peak. How- 
ever, for wavelengths much larger than the typical inter- 
domain distance, a definition of an effective B may be 
consistent with observations. Appendix A shows that if 
only surface viscosities 77^^' and C'^' are important, A" 
has no predicted peak. 

The attenuation effects of elastic and inelastic scatter- 



i ng a re measured by the scattering lengths (|5.8| ), (p.S|), 



(BTQ), and (^.11), as well as the viscous dissipation coeffi- 



cients A" and Im{ki}. Detection of interfacial scattering 
via damping is feasible when the relevant nE is compa- 
rable to A" or Im{ki}. In the transverse — > transverse 
wave scattering, this condition leads to. 



SB_ 
B^ 



> 



P 



,7/41/4 



Tr'^{pa)'^^/i' 



(6.2) 



For a; = 27rx 1000 s"^^ ^ 70 dyne/cm, ly = 0.01 cmVs, 
and a ~ 1 mm, Eq. (Uj) yields n{5BlBof > 1.4. Note 
that this experimentally realizable criterion for detecting 
scattering through wave attenuation depends strongly 
on pa, elastic scattering having a larger effect when the 
wavelength is smaller than the domain radius a. For the 
process involving longitudinal wave scattering pictured 
in Fig. 5(b), the attenuation is proportional to nS/^j. 
The criterion for observing inelastic scattering assisted 
damping is. 



5B 



> 



P 



,5/8^7/2 



i/ii/4CT^'/«7r2(pa)4 



(6.3) 



Though it is p ossible to find parameters such that Eqs. 
( |6.2|) and (6^) are satisfied, the accuracy of the Born ap- 
proximation is questionable when pa or {SB / Bq) are too 
large. However, experiments do show enhanced damp- 
ing when the transverse capillary wavelength is small 
compared to the typical domain size, {pa > 1), and no 
observable d|amping above viscous dissipation for large 
wavelengths .E3 

Alternatively, the scattered waves may be directly de- 
tected. The incident waves vectors are scattered in mag- 
nitude (in the inelastic case) and direction. Using light 
scattering, this experiment may be realized by rotating 
the light scattering plane an angle 9 away from the di- 
rection of the incident wavevector as shown in Figure 6. 
For a small number of elastic scatterers {i/ = 0), the sig- 
nal would be proportional to 1/ (6') P where the scattering 
function is given by Eq. ( |3.6| ). At higher densities, the 
random scattering is expected to effectively contribute to 
the thermal ripplons near the source. A difference in am- 
plitude of the thermal spectra peak measured at 9 ^ tt/2 
is expected depending on whether the wave maker is on. 
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APPENDIX A: INTRINSIC SURFACE STRESS 
EFFECTS 

When surface viscosities are included in the dissipa- 
tive part of the two dimensional stress tensor 11 ^^^^ , the 
behavior of the surface modes can change. For a two 
dimensional Newtonian liquid. 



where r}^"^^ and (^'^^ are the two dimensio nal s hear and 
dilatational viscosities respectively. Using ( |Al[ ) in (p~^), 
the tangential stress balance becomes, 

v±<j{f) + T]^^'Wlvi_ + c*^^Vj_(v_i- t;_L) 

-iy{d,vi_ + V^v,) = f^{r). (A2) 

In Section IV, we neglected surface viscosities. If how- 
ever , B/ lo <^ V^^\C^^\ then the transverse pole is given 
by (EttI) with 



A' 



V2 fcu^y^' x^/'y{y + A) 
6 \(Ja J y'^x + y\f2x + 1 



(A3) 



and, 



A" 



vu: \ y^ \/x/2 + y + 4 



Co / y'^x + j/v 2a; + 1 



(A4) 



(r;( 



2) 



In the above expression, x <^ 1, and y 
(^(2))ti;2/3^j^^i/3 Both A' and A" are monotonic in y. 
The longitudinal mode for uj 3> u'^ /{rj'^'^^+Q^'^^Y becomes, 



h 



{ujv) 



1/4 3iri/8 



and a "shear" mode appears at. 



(A5) 




(A6) 



where the last approximation assumes lo <C v^ jrp'''^ . 

We can also consider the case of a purely elastic inter- 
face; a relevant example may be a polymerized membrane 
on the air/water interface. For an isotropic elastic mem- 
brane, V^ • n (2) is replaced with —5F/5u with the free 
energy of the elastic sheet. 



F==^l d\ [\{r)ul + 2^l{r)u■. 



(A7) 



where Uij = ^{diUj +djUi) is a two-dimensional strain 
tensor and Ui are the in-plane displacements from equi- 
librium. A and /i are spatially varying Lame coefficients 
which may be Ao,/io outside the domains and Ai,/ii in- 
side. With this model, t he tr a nsvers e and longit udina l 
wavevectors are given by (fl.l7|) , (|4.18|) , ([4.19[) , and ( |4.20| ) 
respectively, but with an effective interfacial modulus 
given by 



da 

Bo — 2^0 + Ao — Fo-— 

dl 



r=ro 



(A8) 



n 



'(2) 



r]^'^\d,Vj 



djVi 



SijdkVk) + C, 5ijdkVk, (Al) For /i 7^ 0, an elastic shear mode also occurs at. 




2v^uj 



(A9) 

where the last approxi mati on is va Ud whe n u j ^ u } jv^ . 
The modes given by (pTJ) , (|Il|), (|I|), ( ^!2o|) , and 
( [A9| ) describe a sohd with vanishing shear modulus or 
a liquid membrane with small surface viscosities. Sur- 
face viscoelasticity of the interface can be incorporated 
by allowing B to be complex: B ^ B + iujC. 



ko=P 



c^o 



1/3 



V^¥+^, 



(B4) 



The attenuation of a one dimension al p lane wave is found 
by integrating the y coordinate in (B3), 



{G{r))dy oc e*'=°^ 



where the damping factor a is, 



(B5) 



APPENDIX B: DISORDER AVERAGED 
GREEN'S FUNCTION (G> 

In this section, we develop another means of calculat- 
ing the attenuation of the coherent velocity potential in 
the viscid model. In the limit of small variations in the 
surface tension and sparse uncorrelated scatterers, the re- 
sult is identical wit h th at obtained from the first Born it- 
eration. Equation ( ^.9[ ) can be represented diagrammati- 
cally as in Figure 7(a). The result of ensemble averaging 
over uncorrelated domains, {5(7{r)Sa{r')) ~ eS{r — r'), is 
shown in Figure 7(b). 

In this summation, the cross terms representing corre- 
lated scattering are neglected. Inverting the series and 
defining an ensemble averaged Green's function. 



(G(fc))^Go(fc)£ 



n=0 . 



eGo(<z)fc 



3 / " y 2 



(2^) 



q'\q\Goiq) 



(Bl) 



The cut-off A^^ gives a measure of the typical spac- 
ing between domains of different surface tension. This 
equation is an expansion in the small quantity e. After 
performing the integral over Go (q) and formally summing 
the series, the propagator can written in the form. 



(G(fc))-i = 



CTO 



k\e 



p 



(B2) 



The form of {G{k))~^ is that of Go(fc)"\ but with ctq 
replaced by an effective complex surface tension where 
a I — ep^/6cro and Cfl is a A dependent real part propor- 
tional to e. The Fourier transform of G(fc) yields. 



(G(fc)) 



-le 



2ie 



(ip[<^oi<^R + <^'j)] 



1/3 



H^'\kor) 



(B3) 



where 6 = -^tan 



^(o-j/cTfl) and, 



a — p 



c^o 



.V^'f + '^R. 



1/3 

sin 



-tan ^{ai/aR) 






+ 0{e^) (B6) 



The correlation e for randomly distributed, noninter- 
acting domains can be found by considering the averages. 



and. 



(air))=aoc + ai{l~c) (B7) 



(air)') = ale + a?(l - c) - (a(f))2, (B8) 



where c is the area fraction occupied by the domains. 
Hence, 

{6<j{r)') EE {{a{r) - {air)))') = c(l - c){ao a,)' (B9) 

If we define the correlation on a lattice with spacing of 
the domain size na^ as ((5cr(r)(5cr(r')) = eS{r — f') ~ 
eSij/na^, we find, 



7ra^c(l — c)(o'o — o'l)'^ 



(BIO) 



Using this result in (B6), we find that the attenuation 
a for small area fractions c compares with that found in 
Section III, equations (3.8) and (3.9). 



APPENDIX C: ASYMPTOTIC FORM OF THE 
UNIFORM SURFACE GREEN'S FUNCTION 

We briefly outline the calculation of the asymptotic 
forms of the Green's function matrix for a nonviscous 
compressible two dimensional fluid at the interface of an 
incompressible viscous substrate fluid. The zeroes of the 
determinant of the matrix L((7), 



iB 

(w^ - a\q\q'^ + 2iujvq^){ — q^l + 2vq^ - iu) + 



iB 

q^{aq^ - 2iujvl){ — q^ + 2v\q\) 

UJ 



luq 



(3.1) 



are shown in Fig. 8 and define the modes of the system. 
Note that these eigenmodespare not orthogonal and are 
hence coupled to each otherEll. The Green's function is 
found by Fourier transforming each element of 'L~^{q) via 
the contours shown in Figure 3. Contributions from the 
shear mode near the branch point a.t q = y^iuj/iy and the 
branch cut at are neglected as well as the integrations on 
the imaginary axes. The leading behavior of the Green's 
function is, 



Gii(r) 



G,i(r) 



-UJ 



6pao 



H['\ktr) 



3uji^ en. 



8B^ 



(3.2) 



VLue"^^^^^ (r. 



■AH^\k,r) 
DCTo V r / 

gi7r/8j^5/4^1/4 

^5/2 y n^ 



45, 







{^-^)H^^\hr) (3.3) 



and. 



Gu{r) 



^^TTi/i 



^^('-^)H\^\k,r)- 



6Bo Vr 



4B, 



5/2 



^ Hl^'ikir) (3.4) 



where the coefficients of the Hankel functions are expan- 
sions in small viscosity of the residues at the poles kt and 
ki. Though not needed in this study, similar expressions 
can be found for 6*23 (r) and Gjj{r). 
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FIG. 2. Arrangement for scattering capillary waves from 
circular surfactant domains. The wedge on the left is a wave- 
maker excited at frequency ui. Surface parameters that surfac- 
tant domains change, Xi a-re surface tension a, surface com- 
pressional modulus, B, and interfacial viscosities, r]^^' and 



FIG. 3. A schematic surface pressure - concentration phase 
diagram of an amphiphillic monolayer (cth, is the surface ten- 
sion of a clean air/water interface) . The negative of the slopes 
of the isotherm just outside the plateau multiplied by their 
respective densities, Fj, are the isothermal surface compres- 
sional moduli, Bj , of the two phases near equilibrium. 



FIG. 4. (a) Real and imaginary shifts in the transverse 
wavevector. A' and A" plotted for cto = 72 dyne/cm, u — 0.01 
cm^/s, and uj ~ 5184 s^^. (b) Relative damping rates between 
the transverse and longitudinal modes; the transverse mode 
is more heavily damped for parameters above the y = axis, 
longitudinal waves are damped quickly below y = 0. 
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FIG. 5. (a) Limit when longitudinal mode is damped FIG. 7. (a) Diagram representing the sum in (Al). (b) 

near each scatterer. The transverse waves are weakly Green's function averaged over frozen disorder, 
damped on interdomain length scales I. (b) Longitudinal 
waves are long-ranged, transverse waves are strongly damped; 

Im{kt} <C l^^ ■ FIG. 8. Approximate location of poles in the complex 

g-plane. The shear mode in our analyses has a large imaginary 
component. The contours are the integrals used to evaluate 

FIG. 6. Schematic of a proposed experiment to detect scat- the Green's functions 
tering from interfacial inhomogeneities at the air /water inter- 
face 



11 



